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Quantum Lyapunov control was developed in order to transform a quantum system from arbitrary 
initial states to a target state. The idea is to find control fields that steer the Lyapunov function 
to zero as t — > oo, meanwhile the quantum system is driven to the target state, fn order to shorten 
the time required to reach the target state, we propose two designs to optimize Lyapunov control 
in this paper. The first design makes the Lyapunov function decrease as fast as possible with a 
constraint on the total power of control fields, and the second design has the same purpose but with 
a constraint on each control field. Examples of a three-level system demonstrate that the evolution 
time for Lyapunov control can be significantly shortened, especially when high control fidelity is 
required. Besides, this optimal Lyapunov-based quantum control is robust against uncertainties in 
the free Hamiltonian and decoherence in the system compared to conventional Lyapunov control. 
We apply our optimal design to cool a nanomechanical resonator, a shorter cooling 
time is found with respect to the cooling time by the conventional Lyapunov design. 

PACS numbers: 03.65.Yz, 02.30.Yy, 03.67.-a 



I. INTRODUCTION 

Quantum control [jj, y] has attracted much attention 
in recent years and it has found potential applications 
in many fields such as quantum information processing, 
quantum chemistry and quantum simulation. Among 
these quantum control problems, state transfer is a cen- 
tral task. Various methods such as quantum optimal 
control and Lie group decompositions have been used to 
design control laws to drive quantum systems to target 
states or to realize some specific operations [3h13||. 

Quantum Lyapunov control was proposed in the early 
2000s as a good candidate for state transfer [Tj,[l3- This 
strategy has been widely studied recently both in the- 
ory and applications 16-26], because it offers a simple 
and effective way to design control fields. However, the 
problem of speeding up Lyapunov control has not been 
widely considered to date. Quantum Lyapunov control is 
practically employed in an open-loop way without mea- 
surement and feedback. Hence, it is quite natural to 
shorten the evolution time so as to overcome the deco- 
herence effect induced by inevitable interactions with the 
surrounding environment. 

In quantum Lyapunov control, a function V, called a 
Lyapunov function of quantum states, is specified to de- 
sign time- varying control fields. The system converges to 
the target state given by V — while V decreases to its 
minimum. Based on this concept, we present a scheme to 
optimize Lyapunov control by using the following idea: 
Speed up evolution to the target state by making V de- 
crease faster. 

Generally, the total Hamiltonian in a quantum control 
system can be written in two parts Hq + Yl n =i fn(t)H n , 



where Hq is the free (internal) Hamiltonian and H n are 
external control Hamiltonians with /„ (t) representing the 
corresponding control fields. These control fields should 
be designed to ensure V < 0. In fact, there are many 
ways to choose f n (t) to achieve this goal. In this paper, 
we consider the question: With given control Hamiltoni- 
ans H n , how should we design the shape of the control 
fields f n (t) m order to make V decrease fastest. 

It is shown that one can enable V to decrease faster 
simply by enhancing the strength of the control fields. 
However, strong control fields may not be feasible and 
always bring unwanted results, for example a large en- 
ergy (power) cost, invalidation of mathematical approx- 
imation and treatment for the system. In view of these 
factors, we propose two designs for the control fields. 
One is under the constraint that the power-type quan- 
tity W — J2 n =i fn{t) 2 is bounded. The other is under 
the constraint that the strength of each field is bounded. 
These designs for control fields make V decrease as fast as 
possible within given limitations. The second design has 
the simple form of "bang-bang" control which is easy to 
be implemented in an experiment. We also illustrate our 
control method with a three-level system. The results 
suggest that the evolution time is significantly reduced 
compared with the conventional method. 

The paper is organized as follows: In Sec. II, we de- 
rive control fields for a general Lyapunov function V — 
Tr(Pp) under the two aforementioned constraints. In 
Sec. Ill, we simulate our field designs and compare them 
with the conventional method. The robustness of the op- 
timal Lyapunov designs is analyzed in Sec. IV. Finally, we 
summarize our work in Sec. V. 



II. DESIGN OF CONTROL FIELDS 

In quantum Lyapunov control, the system is steered 
from an initial state to a target state by control fields 
designed using Lyapunov function V. The goal of this 
paper is to obtain optimized control fields that make the 
time derivative of the Lyapunov function \V\ largest so 
as to speed up the evolution to the target state. We start 
from a closed quantum system described by the Liouville 
equation 



dp 
~dt 



-i[H + H c (t),p], 



(1) 



where Hq is the free Hamiltonian for the controlled sys- 
tem and H c (t) is a time-dependent Hamiltonian repre- 
senting coupling to external control fields which is called 
the control Hamiltonian. We have set h = 1 and assume 
that the system is controllable. In Lyapunov control, the 
solution of Eq.([T]) converges to the minimum of V{p). 
Meanwhile, the state converges to a set of states char- 
acterized by the La Salle's invariance principle 1]. The 
control Hamiltonian H c (t) can be written in the form, 



H c {t)^Y,fn{t)H n 



(2) 



n=l 



where H n (n = l,...,fc) are time-independent Hcrmitian 
operators corresponding to different types of external 
control and /«(£) are time- varying real functions, usu- 
ally representing electro-magnetic fields, k is a positive 
integer. 

In this paper we consider the following form of Lya- 
punov function 



V = Tr(Pp), 



(3) 



where P is a Hermitian operator and assumed to be posi- 
tive semi-definite in order to satisfy the standard require- 
ment for a Lyapunov function, V > |2|. Also, some 
other forms of Lyapunov function can be described by 
Eq.([3]), such as that based on the Hilbert Schmidt dis- 
tance diii]. 

The time derivative of V needs to be calculated to 
design the control fields, 



V = Tr(-iP[H + ]T f n {t)H n ,p\) 



n=l 

k 

Tr(-ip[P, H ]) + Y, fn(t)Tr(-ip[P, H n \) (4) 



k 

E 

n=l 



fn(t)T n 



where T n = Tr(— ip[P, H n ]) is a real function of p, H n and 
P. We have used the assumption that [P, Ho] = 0, which 
can be achieved by constructing P using the eigenvectors 
of -ff . 



The Lyapunov control strategy requires V < 0. There 
are many ways to design f n (t) to satisfy this requirement. 
A simple and conventional way is to let f n (t) = —KT n 
with K > so that 



y(p) = -]T^T n (i) 2 <0. 



(5) 



With such control fields, V will decrease to its minimum 
and state p will converge to the target state pf with 
the same spectrum as the initial state po and satisfying 
Tr(e- iHot p f e iBot [P,H n ])=0 0. 

In the conventional field design method, the ampli- 
tude of the control fields /„(£) is proportional to T n . 
That means when T n is small (for example, when p is 
very close to pf), f n (t) will become small leading to a 
slow decreasing of V and a long evolution time. Our aim 
is to determine optimized control fields f n (t) to enable 
V to decrease as fast as possible in order to speed up 
control. From Eq.Q it is seen that if each f n (t) has a 
different sign to the corresponding T n , then V < and 
large |/„(t)| will lead to fast decreasing of V. Therefore, 
the problem has to be discussed under a constraint on 
the control fields f n (t). Considering the following rea- 
sons for constraining the control fields: First, one often 
wishes control fields to be weak in order to reduce energy 
(power) costs. Second, strong external control fields may 
lead to invalidation of the modeling of the system. Third, 
strong fields may disturb neighboring quantum systems 
that we do not want to disturb. We propose two designs 
of control fields under constraints on the power of the 
control fields (constraint A) and on the strength of each 
control field (constraint B), respectively. 



A. A power-type Constraint 

First, we consider the following constraint on the con- 
trol fields 



k 



(6) 



Since the control fields /„(£) are always associated with 
the amplitude of the electro-magnetic fields, the quantity 
W can be interpreted as a power-type quantity. We will 
call it power for simplicity in the following. The total 
power of the control fields is bounded in this case. 

Consider J2 n=1 T% ^ (p / pf) in time t and the 

constraint for f n (t) is J2 n =i /«(^) 2 = W. In order 
to determine the optimized control fields that minimize 
V = X) n =i fn(t)T n (V is negative), we use the Lagrange 
multiplier method. Let 



L = J2fnTn+Hj2^~W) 



(7) 



where A is the Lagrange multiplier, and /„ represents 
f n {t) at a certain time t. Then from the following equa- 
tions 



■^-L = T n + 2A/„ = 0, 

dfn 



(8) 
(9) 



it is easy to obtain the amplitude of the fields f n at time 
t, 



In ■ 



WT n 






(10) 



with T n = Tr(— ip[P, H n \). The corresponding time 
derivative of Lyapunov function is 



v, = 5>T„=-^£kS 



n=l 



(11) 



EK rp2 

n=l n 



It is seen that V is proportional to VW, so we choose 
W = W max for a faster decreasing of V and our control 
design for all evolution time reads 



•JWTr 



fn(t) 



Vzi 



T'i 



(ELl^nVO) 



o 






0). 



(12) 



Note that when p reaches the final state pf, all T n = 
Tr(—ipf[P,H n ]) become zero and all control fields are 
switched off. Considering that p converges to pf asymp- 
totically, we will switch off the control fields after 
D(p, pf) < e where D denotes some measurement for the 
distance between p and pf and e is the required precision. 
In the case k = 1, i.e., there is only one control Hamil- 
tonian, our control design reduces to 



hit) 



-Vw„ 



o 



(Ti > 0) 
(Ti < 0) 
(Ti = 0) 



(13) 



which has a simple "bang-bang" control form. With its 
discrete shape of the control fields, this control design 
should be easy to realize experimentally [g, [23| ■ 



B. Constraint on the strength of each control field 

Next, we will find the optimized control fields when 
the strength of each field is bounded. For simplicity, we 
assume the maximum strength of every control field /„ (t) 
is S (S > 0), i.e., 



\fn(t)\<S, (n= 1,2, •••,£). 



(14) 



From V = Yl n =i fn(t)T n , it is easy to obtain the op- 
timized control fields that minimize V with condition 



Eq.jH 



/»(*) = 



■s 


(T„ > 0) 


s 


(T n < 0) 





{T n = 0) 



n=l,2,---,k (15) 



and the time derivative of the Lyapunov function is 



V=J2fn(t)T n = SY,\Tn\- (16) 



This design has the "bang-bang" control form with k 
different control fields. When there is only one control 
Hamiltonian, the design has the same form as that in 
Eq.CdD with VW replaced by S. 

We have presented two designs of control fields for sys- 
tems described by Eq.([T]) and Lyapunov function Eq.Q 
with two constraints. However, these designs can also be 
applied to other Lyapunov control as long as the deriva- 
tive of Lyapunov function has the form of Eq. Q where 
f n (t) represents a control field and T n represents a real 
function of the quantum state. In fact, for many differ- 
ent kinds of Lyapunov function and different dynamical 
equations [3 |jjt HO, HJ , V takes this form. 



III. ILLUSTRATIONS 

In this section, we will present an example to illus- 
trate the proposed schemes. The example consists of a 
3-level system driven by a control Hamiltonian. We show 
that the system can be steered to an eigenstate of the 
free Hamiltonian from arbitrary initial states (except the 
states in the La Salle's invariant space) by both the con- 
ventional design f n {t) = —KT n and the design proposed 
in this paper. The difference is that the present design 
can speed up the convergence. 

Consider a three-level system described by the quan- 
tum Liouville equation 



^=-i[#o + ^ /„(*)#„, p] 



with free Hamiltonian 



Hn 




(17) 



(18) 



where the energy difference between |1) and |2) 

(|2) and |3)) is cj (^w). H = 1 has been set through- 
out this paper. 

The aim is to steer the system from an arbitrary initial 
pure state po = \(po)(4>o\ to an eigenstate, say \4>f) — 
1 3) = [1,0, 0] T of the free Hamiltonian H$. We choose 
the Lyapunov function V = Tr(Pp) with 



P 




(19) 
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FIG. 1: (color online) (a) Pf (solid line) and \V\ (dashed 
line) as a function of time with the conventional control de- 
sign f n (t) = —KT n . (b) Control fields (in units of to) for 
W mO B = 0.0001, and S ma x — 0.007 . These fields are switched 
off when P } > 0.999 at t = 191 . 
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FIG. 2: (color online) (a) P f (solid line) and \V\ (dashed 

line) as a function of time with the control design Eq. (|12[) 
under constraint A. (b) Control fields (in units of w) with 
Wmax — 0.0001. These fields are switched off at t = 90 when 
P f > 0.999. 



According to the Lyapuonv control theory, the system 
Eq.(fTT|) will be driven to an eigenstate of P with the 
minimum eigenvalue 0, i.e., p — > pf = |3)(3|. Addi- 
tionally, in the pure state case, this Lyapunov function 
can be explained as the Hilbert Schmidt distance be- 
tween \4>) and |0/). Recall that P = I — \<f>f)(4>f\, and 
V = Tr(Pp) = Tr((I - |0/)(0/|)|0)(0|) = 1 - \(4>\4>f)\\ 
our control design in this example acts to make the dis- 
tance between \<f>) and \<j>f) decrease as fast as possible 
with a given restriction on the control fields. 

In order to achieve the best performance, we choose 
the control Hamiltonians H n to be the generators of 
the SU{2>) group A„(n = 1,...,8), namely, H c (t) is con- 
strained from 



Ai = 



control Hamiltonians are chosen as, 



A, = 




(20) 



A 7 = 



Notice that since only Ai, A2, A4 and A5 satisfy [P, A„] 7^ 
0, only these generators are effective for our model which 
can be understood by examining Eq.(j3]). Therefore, the 



H 1 = X 1 = 



Ha 




• Ho — At — 



, Ha — Ah 




which can be rewritten as H\ = |3)(2| + |2)(3|, H2 = 
t(-|3)<2| + |2)(3|), H 3 = |3)(1| + |1)(3|, and H A = 
i{— 13)(1| + |1)(3|), which couple the energy levels |1) 
and |2) to the final state |3). We would like to note 
that the control Hamiltonians in this example are opti- 
mal, because any operator for this 3-level system can be 
written as an expansion of these generators. For high- 
dimensional system, the problem becomes more compli- 
cated. 

We first simulate the problem with the conventional 
control field design f n (t) — —KTr(—ip[P,H n ]) with 
K = 0.01 and initial state \<j> ) = A=(\l) + |2) + |3)). 
The control fields are switched off when the probability 
Pf = \(<j)\<j>f)\ 2 reaches 0.999 both in this simulation and 
the following two so as to compare the evolution time. 
FigOJa) shows the evolution of probability Pf (black 
solid line) and \V\ (blue dashed line). It is seen that the 
system is driven to target state |</>/) and the evolution 
time is about t = 191 for Pf = 0.999. The time-varying 
control fields f n (t) are plotted in Figfljb). The power 
W = Yl n =i /«(*) 2 reaches its maximum W max = 0.0001 
at t = 17 and the maximal strength of a single control 
field is |/i I = 0.007 at t= 16. 
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FIG. 3: (color online) (a) Pf (solid line) and \V\ (dashed 

line) as a function of time with "bang-bang" control Eq. (|15[) 
under constraint B. (b) Control fields (in units of u> ) with 
S = 0.007. These fields are switched off at t = 74 when 
P f > 0.999. 
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Next, we employ the control design Eq. (jT2|) under the 
constraint on the power of control fields. In order to be 
comparable with the last example, the maximal power 
W max and initial state are chosen to be the same as the 
above one. Results in FigJ^fa) show that the evolution 
time for Pf = 0.999 is t = 92 which is evidently shorter 
than that shown in Fig. [TJ Control fields are plotted in 
Fig|2jb) which have a sinusoidal shape. In this example, 
for all the time before control fields are switched off, the 
power of control fields remains constant and the shapes 
of fields are optimized so that the Lyapunov function 
decreases fastest under constraint A. We can see this from 
evolution of \V\ in Fig|2fa) (blue dotted line). 

Now we study the control design Eq. (fT5"j) under con- 
straint B. For comparison, we let the maximal strength 
S and the initial state be the same as that in the first sim- 
ulation. The simulation results are illustrated in FigJ3] 
The evolution time for Pf = 0.999 is t = 73 which is 
about 38 percent of the first example. In this exam- 
ple, the control fields shown in Figf^b) are step-like. 
Such a control method makes Lyapunov function de- 
crease fastest under the restriction of strength and has 
the advantage of being easy to be implemented in ex- 
periment. The evolution of \V\ is not smooth (shown in 
FigHUa) by blue dashed line), which is due to the discrete 
control fields. 

Furthermore, we plot in FigU the evolution of 
D = 1 — | (</>|</>/) | 2 for the three designs of control 
field with 50 randomly chosen initial states ( \<j) r ) — 
R[ne l27Tri ,r 2 e l27Tr!i , r 3 e l27rr6 ] T , where n(i = 1,...,6) are 
random numbers uniformly created between and 1, 



FIG. 4: (color online) Evolution of D = 1 - \((f>\</>f)\ 2 for the 
three control designs with logarithmic scale. (a),(b) and (c) 
correspond to the conventional control field design, the design 
under constraint A, and that under constraint B, respectively. 
Each picture is a result averaged over 50 random initial states. 
The convergence speed of (b) and (c) is faster than that of 
(a) as these figures show. 



R = ===== is a normalization factor) respectively. 

yrf+ri+ri ; P J 

The simulations show that the convergence rate for the 
usual control design is exponential where the evolution 
time grows linearly with the distance between the actual 
and target states D, whereas the convergence rate of our 
two methods is larger than the conventional one, espe- 
cially when the control fidelity is required to be high. 
The reason is our methods keep \V\ at its maximum for 
all the evolution even if T n is very small. 

We found that some of the control fields may oscil- 
late with very high frequency at the end of the control 
in Ea. (fT5|) with constraint B. The oscillation depends on 
the initial state and the distance between the actual and 
target states. The reason of this oscillation is as fol- 
lows. When the state is close to the target state, some 
T„ become very small leading to ineffectiveness of the 
corresponding control field f n {t). While in conventional 
Lyapunov control, the control field /«(£) decreases to zero 
with T n — > 0. In our method, however, f n (t) is designed 
to take the value S or —5, which makes the state oscil- 
late almost every step of simulation. This problem can 
be solved by averaging the control fields over a proper 
time period and then use the reshaped fields instead of 



the oscillating one. In fact, this average can be used in 
the situation when the control fields are not a (fast) os- 
cillating function of time, it yields the same control fields 
and maintains the results. 

It is worth noting that the convergency depends on the 
amplitude of the control fields, when the control fields os- 
cillate very fast, the Lyapunov function V almost stops 
to decrease. This is different from the conventional Lya- 
punov design where V keeps decreasing with the ampli- 
tude of the control fields approaching zero. This suggests 
that we can adjust the convergency by properly designing 
the amplitude of the controls fields. We will show and 
prove this point in the future work. 



IV. ROBUSTNESS OF THE DESIGNS 

One may wonder if these optimal designs improve the 
robustness of the Lyapunov-based control. In the follow- 
ing, we shall examine this problem following the repre- 
sentations in [26J by calculating an average fidelity of the 
system with Hamiltonian uncertainties, decoherence and 
field fluctuation in the controls Eq.lfTT]). Here the aver- 
age fidelity is defined by P' f = ^(E^Li \{^M})?)^ \<t>j) 
denotes the actual state evolving from a random initial 
state under the control with uncertainties, decoherence 
or fluctuations. In other words, the average is taken over 
N actual states, each evolves from a randomly chosen 
initial state, driving by the controls with uncertainties, 
decoherence or fluctuations. Our focus is on whether 
the optimal designs is robust against these uncertainties 
compared with the conventional one. 

We begin with analyzing robustness against the uncer- 
tainty in the free Hamiltonian Hq. The uncertainties can 
be taken into account by adding a perturbation SHq to 
the free Hamiltonian, i.e., 

H o — >• H + SH . 

Here SHq — E n =i &n^n with 5 n a real number and \ n 
the generators in Eq. (|20[) . For simplicity, we examine 
separately the 8 uncertainties SHq = S n X n , n = 1, 2 • • • 8. 
The equation f t = -i[H + SH + £* =1 f n (t)H nt p] is 
simulated for the three designs with the same control 
fields and the same parameters as in section III. Selected 
results are showed in FigJS] where the fidelity P f is an 
average over fidelities from 1000 randomly chosen initial 
states. 

The simulations show that, (1) the optimal Lyapunov 
control is robust against the uncertainty Ai, and is sen- 
sitive to that of A3; (2) the bang-bang design is more 
robust than the design with power constraint, and the 
conventional design has the worst robustness. In fact, 
the robustness against the uncertainty <5gAs is simi- 
lar to FigfSJa) and robustness against the uncertainties 
S,An(n = 2,3,4,5,6,7) is roughly similar to Fig[5jb). 
Thus here we only show the robustness of the control 
against the uncertainties Ai and A3. This conclusion con- 
firms that the system is sensitive to the uncertainties that 
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FIG. 5: (color online) The average fidelity verses uncertain- 
ties in the free Hamiltonian, (a) for Ai and (b) for A3. The 
blue triangle, green star and black circle represent the fidelity 
obtained by bang-bang control design, design with power con- 
straint and the conventional one, respectively. 
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FIG. 6: (color online) The average fidelity as a function of 
the decoherence rates, here we choose 72 — 73 — 7. The blue 
triangle, green star and black circle represent bang-bang con- 
trol design, design with power constraint and the conventional 
one, respectively. 



commute with the free Hamiltonian, while it is robust 
against the uncertainties that does not. 

We next discuss the robustness against decoherence. 
In this discussion, we assume that states |2) and |1) are 
long lived, namely only the spontaneous emission |3) — > 
|1) and |3) — > |2) is assumed. This can be described 
by C{p) = E^aTifo" P°t ~ \ G ^P°t - l^P^t) with 
G I = N)(3| an d o1[ — |3)(i|. With this assumption, 
the target state is not a steady state, so the fidelity will 



be affected seriously. However, the optimal controls are 
better than the conventional one due to the short time 
needed to drive the quantum system from an initial state 
to the target state. As we did in FigO we calculate the 
average fidelity for the three control designs with 1000 
random initial state and fixed 7. The results are depicted 
in FigJSJ In this case, the fidelity of the two optimal 
designs is obviously higher than the conventional design. 
The robustness against field fluctuations (with zero- 
mean) and errors in the initial state is also explored. The 
results are similar to that in [26| , namely the control field 
fluctuation with zero mean affects the fidelity slightly, 
while it depends sharply on the fluctuation with non-zero 
mean, and the final fidelity is sensitive to the initial state. 
In these cases, the optimal control has no advantage with 
respect to the conventional one. 



APPLICATION TO THE COOLING OF A 
MECHANICAL OSCILLATOR 



The cooling of mechanical resonators [27H31J becomes 
an active research topic in recent years due to its poten- 
tial applications in detecting extremely small displace- 
ment and observing quantum phenomenon of macro- 
scopic mechanical object. In this section, we apply 
the Lyapunov control to cooling a nano-mechanical res- 
onator. The results show that it is possible to cool a 
nano-mechanical system to its ground state by Lyapunov 
control, the optimized control design leads to a shorter 
cooling time with respect to the conventional control de- 
sign. 

Consider a nano-mechanical resonator (called target) 
with frequency uj coupled to the other microwave (opti- 
cal) oscillator (auxiliary system), the microwave oscilla- 
tor has a sufficiently higher frequency J7 such that it can 
be prepared in its ground state at finite temperature. In 
the language of Lyapunov control, the free Hamiltonian 
of the composite system is given by 



Ho = Huja [ a + Mlb'b. 



(22) 



We assume the coupling Hamiltonian of the two oscilla- 
tors has the following form, 



H c = g(t)xAXB 



(23) 



with xa = a + a) and ig = b + tf . This type of Hamil- 
tonian can be realized by coupling the target to a LC 
oscillator and the coupling rate g(t) can be modulated 
by the voltage of the LC circuit |27h29| . 

In the sideband cooling, g(t) is modulated at $7 — uj so 
that the two resonators are effectively coupled and the 
rotating wave approximation (RWA) applies when the 
coupling g is weak. Recently, the authors of 27] shown 
that quantum control can improve the cooling, when the 
control goes beyond the RWA in the ultra-strong coupling 
regime g ~ u> [23] ■ Here, we show that we can obtain 
the control design by the Lyapunov functional, and the 
optimized Lyapunov design can shorten the cooling time. 
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FIG. 7: (color online) (a) Phonon number versus time for the 
conventional Lyapunov design (upper half) and optimal Lya- 
punov design (lower half). The evolution time for (n a ) = 0.05 
is tot — 37.5 (upper) and cut — 13.4 (lower), respectively, (b) 
Time-dependent coupling for the conventional design (upper) 
and the optimal design (lower), where g(t) for the two designs 
share the same maximal strength \g m ax\/uj = 1.91. g(t) for 
the optimal design is shut off after cot — 13.4. 



Denote the state of the two resonators by p, we can 
choose the Lyapunov function as 



V(p) = Tr(a f ap) = 



(24) 



namely, we choose the mean phonon number of the tar- 
get resonator as the Lyapunov function, which is non- 
negative and becomes zero when the target system is 
cooled to its ground state. By the same procedure, we 
get 



V(p) = g(t)T n . 



(25) 



Here T n = Tr{-ip[rfa, (a+ + a){tf + b)}). If g(t) is 
chosen to keep V < 0, the phonon number of tar- 
get resonator will decrease monotonically. The conven- 
tional Lyapunov design for the time-dependent coupling 
strength is g(t) = —KT n with K a positive constant, 
while the optimal Lyapunov design is a bang-bang con- 
trol given by Eq. (fT5|) . 

Assume that the target resonator is initially in a ther- 
mal state with average phonon number (n a ) =6.38 and 
the auxiliary system is prepared in its ground state. Note 
that the frequency of auxiliary system $7 should be suffi- 
ciently large compared with uj, such that thermal fluctu- 
ation has small effect on the microwave oscillator. Our 



numerical simulations show that larger f2 leads to faster 
oscillation of g(i), but it does not affect the cooling re- 
sults. Here we set O/w = 20. 

The simulations are performed in the Fock space, we 
truncate the Fock space of each oscillator up to 20 Fock 
states and the dissipation of each resonator is ignored. 
Compare to simulations with 25-Fock state-truncations, 
improvement is not significant, so the simulations with 
20-Fock state-truncation are reasonable. We compare 
the evolution of phonon number (n a ) and g(t) for the two 
designs in FigfTJ The top half of FigJ7][a) represents the 
evolution of (n a ) for conventional design with K = 0.03 
and maximal control field strength \g(t)\/oj = 0.191. It 
takes t — 37.5 (in unit of 1/w) for the target to reach 
(n a ) = 0.05. In contrast, for the optimal design, the 
evolution time for reaching the same phonon number is 
t = 13.4 as shown in the lower half of Fig|7]Ja). Obviously 
the optimal Lyapunov design shorten the cooling time. 
FigE^b) shows the time dependence of g(t) for the con- 
ventional (upper half) and optimal (lower half) design. 
For each case, the control field g(t) starts with a non-zero 
small number to avoid T n = at the initial time. It is 
seen that the optimal Lyapunov design leads to a faster 
decrease of Lyapunov function (i.e. the phonon number) . 
In addition, the major components of the oscillation fre- 
quency of g(t) is automatically turned to the frequency 
difference O — uj of the two resonators (like that in the 
sideband cooling scheme). 



VI. SUMMARY 

We have presented two designs for Lyapunov con- 
trol under constraints on the total power and individ- 
ual strength of the control fields. These designs make 
Lyapunov function decrease fastest determined by the 
constraints. It has been shown that the implementation 
of our designs leads to a shorter time towards the tar- 
get state especially for high fidelity requirement. More- 
over, the second control design gives simple bang-bang 
control fields, which may be easy to implement in exper- 
iment. Intuitively, our methods use a constant power or 
strength of control fields to make Lyapunov function de- 
crease as fast as possible. This optimal control is more 
robust against uncertainties in the Hamiltonian and de- 
coherence in the system with respect to the conventional 
design. We also explore the application of our op- 
timal design to cool a nanomechanical system, a 
significantly shorter time is obtained compared 
with the conventional Lyapunov design. Here we 
focused only on how to design control fields with fixed 
control Hamiltonian H n , a general formalism to choose 
control Hamiltonian H„ to speed up Lyapunov control is 
still an open issue. 
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